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The universal behavior of the continuous superconducting (glass) phase transition is studied in 
a YBa2Cu307_a thin film. A new analysis technique for extracting critical exponents is developed, 
using current-voltage (J-E) measurements. The method places narrow limits on the glass scaling 
exponents, z g and v g , as well as the glass transition temperature, T g (H). The universal values 
z g = 3.8 ± 0.2 and v g = 1.8 ± 0.1 are obtained. Using these values, the data are scaled for fields 
in the range 1-9 T. Collapse of the scaling forms at different fields describes the universal scaling 
function. The field dependence of the universal collapse is interpreted in terms of the 3D XY 
multicritical scaling theory. 
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I. INTRODUCTION 



The nature of the superconducting phase transition at 
finite fields in high temperature superconductors [HTS] 
remains surprisingly elusive. In clean systems, recent 
evidence suggests a first order melting transition of the 
vortex latticeE However in disordered systems, the issue 
is not completely resolved: Does superconducting order 
exist at finite temperatures, separated from the normal 
state by a well-defined continuous phase transition? Or 
do manifestations of superconductivity at nonzero tem- 
peratures represent only nonequilibrium phenomena, as 
in the flux creep descriptionTa This question is of great 
practical interest since the two mutually exclusive de- 
scriptions involve distinct dynamic equations which gov- 
ern the response of driveruvortices. In one case, the equa- 
tions are mean-field-like;u in the other, they describe a 
(glassy) critical state, which can occupy a fiery wide re- 
gion of the superconducting phase diagram.! From a fun- 
damental viewpoint, the debate encompasses the nature 
of the superconducting state and long range order in the 
cuprate superconductors. 

Shortly after the discovery of HTS, significant progress 
was rnacLe to clarify these points. Several types of glassy 
statesErQ were proposed as alternatives to the creep pic- 
ture. A key feature of these glass phases is the true 
absence of resistance at low currents, for temperatures 
below the glass transition line, T g (H). Initial experi- 
ments strongly supported the presence of such a contin- 
uous phase transitionJjEl and the superiority of the glass 
description over the creep picture.£3 However, further 
progress is hampered by the absence of a firm descrip- 
tion of the glassy state, particularly regarding the nature 
of the glass transition. Neither the critical scaling expo- 
nents nor the scaling functions have been determined an- 
alytically, due to the difficulty of the theoretical problem. 
Further experimental treatments, while individually sup- 
porting a scaling descriptionO of current-voltage (J-E) 
characteristics, appear as a group to violate the notion of 
critical universality: both sample- and field-dependence 



in the scaling exponents are observed by some authors, til 
A more rigorous test of universality, involving the scal- 
ing functions (which must also be universal according 
to the critical theory), has been investigated previously 
by Xenikos et raft3 in terms of sample-dependence, by 
Woltgens et a/£3 in terms of the angle dependence of H, 
and by Kotzler et alSj in terms of magnetic field depen- 
dence. Since further theoretical advancements are not 
immediately forthcoming, confirmation of the existence 
of true finite temperature superconductivity now rests 
on the ability to discern universal behavior near the pur- 
ported glass transition. We note, however, that this test 
provides evidence only for critical scaling, and gives lim- 
ited insight into the nature of the glassy state. 

Other types of critical fluctuations (besides glass fluc- 
tuations) may arise and compete with the glass exci- 
tations for critical dominance. If the scaling regions 
overlap, the corresponding scaling theories must be self- 
consistent, a fact which has a profound implications for 
universality. One such competing excitation in HTS is as- 
sociated with the 3D XY model. The critical behavior of 
the zero field superconducting transition, T c , is expected 
to belong to the 3D XY model class, in the strongly type- 
II limit.El Furthermore, the finite field superconducting 
transition line T g (H) terminates at the 3D XY multi- 
critical point: T c = T g (H = 0). The topology of the 
H-T phase diagram therefore ensures that the glass and 
3D XY critical regimes will overlap, and that their uni- 
versal behaviors will be intimately pelated in some region 
about the 3D XY transition point .E£rE3 The cautions ap- 
plied to any mean-field dynamic description in the glass 
scaling region also apply to the 3D XY region, vthjeh. 
extends over a wide temperature range at4ew fields 
encompassing the low field glassy region.ll3iI3 

In this work, a single YBa2Cri307_5 thin film is stud- 
ied. In order to explore universal behavior, a rigorous 
scaling treatment is required to place narrow limits on 
the glass scaling parameters. After restricting the anal- 
ysis in this way, a J-E scaling procedure is performed, 
adjusting the scaling parameters within their bounds to 
obtain the best results. To check universality, glass scal- 
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ing is compared for different fields, not only in terms of 
scaling exponents, but also scaling functions; both are 
found to be field independent. However, fields do enter 
into the glass scaling in a way which is determined by the 
3D XY multicritical theory at low fields. Although the 
present treatment is not applied to different samples here, 
it is hoped that such an extension may be performed in 
the near future. The following analysis is presented in 
the vortex glass formulation,El however this does not con- 
fine the critical analysis in any way, since the different 
J-E scaling formulations are isomorphic to one another. 



II. EXPERIMENT 

A 4000 A thick, deoxygenated YBa 2 Cu 3 07„5 film, 
with the c-axis oriented perpendicular to the film surface, 
was used in this study. The film was prepared by pulsed 
laser deposition onto a heated (001) LaAlC>3 substrate, 
and initially had a superconducting transition tempera- 
ture of T c — 92 ± 2 K at zero applied field, ft was then 
deoxygenated at 400 °C using a controlled flow of oxy- 
gen and argon gasses, followed by a rapid quench to room 
temperature. The resulting superconductor had a transi- 
tion temperature T c — 77 K, corresponding to an oxygen 
stoichiometry of 5 ~ 0.24. The film was patterned into a 
120 /im x 3100 (im bridge, using laser ablation. 

The resistivity, p(T), was measured as a function of 
temperature for fixed currents and fields. Measurements 
were then performed for a range of applied current val- 
ues using a conventional four point geometry, for mag- 
netic fields in the range 1 < H < 9 T. The current 
densities ranged from 20 to 4 x 10 4 A cm~ 2 . The tem- 
perature range over which the resistivity was measured 
depended on the current. Typically, data were collected 
on warming for temperature windows of 25 K, starting 
at temperatures where the voltage signal increased over 
the instrumental noise level. This ensured that data were 
collected on both sides of the glass transition tempera- 
ture, T g . Throughout this work, errors are quoted as 2a 
(standard deviation) unless otherwise noted. 



III. SCALING ANALYSIS 



A. Determination of z g and T g 



The main focus of this work is the universal critical be- 
havior near the glass transition temperature, T g . There- 
fore a systematic, precise method for determining T g , and 
the glass scaling exponents z g and v g is required. One 
method for identifying these parameters involves scal- 
ing of current-voltage (J-E) isotherms, obtained in 
vicinity oLT g for a single applied magnetic field (H) 
TypicallytiJ the three scaling parameters are varied si- 
multaneously to obtain the best collapse of the isotherms 
With no constraints placed on the fitting parameters 



this method can obtain a reasonable scaling collapse for 
a large range of exponents. Conventionally, constraints 
are placed on the fitting procedure by identifying z g and 
T g from the critical isotherm,Q then requiring that the 
final scaling adjustments do not stray far from this es- 
timate. However, low densities of isotherms, difficulties 
in subtracting the normal state background (causing the 
critical isotherm ta appear as an imperfect power law), 
and Joule heatingjj all create ambiguities in the identifi- 
cation of the critical isotherm, reducing the precision of 
T g to 1 K or worse. 

To improve this situation, an alternative scaling proce- 
dure is developed below, in which z g and v g are system- 
atically obtained. This procedure places more rigorous 
constraints on the fitting parameters by largely circum- 
venting the difficulties described above. An important 
feature of the method is its optimization of scaling in the 
vicinity of T = T g \ although this condition is required 
theoretically, it is not rigorously enforced in most scaling 
techniques. A similar technique was applied previously 
in the study of 3D X Y_ low field fluctuations in high-T c 
superconducting films.li3 In the present case, the starting 
point for the scaling procedure is the following ansatzfl 



a(JfT)^- 1 ^ 2 = F \t(J /T)- 1 / 2 "* 



(1) 



Here, a — 1/p = J/E is the nonlinear conductivity, F[x] 
is a scaling function, and t — (T — T g )/T c is the relative 
temperature variable, for a fixed field H . 

Since the left hand side of Eq. (Q) displays no singulari- 
ties, the function F[x] must be smooth and single- valued. 
In particular, it is analytic at the glass transition T = T g , 
where its argument x — t(J "/J 1 )" 1 / 2 " 9 (the vortex glass 
scaling variable) is equal to 0. When t = 0, the scal- 
ing variable no longer depends on the current density 
J; both sides of Eq. (|l|) should then be independent of 
J, provided that z g is chosen correctly. The resistivity 
curves corresponding to different current values therefore 
all cross at T g when plotted as in Fig. 1. In the inset, 
the exponents z g , which provide the best simultaneous 
crossing, are shown for each applied field. The error bars 
reflect the fact that satisfactory crossings are obtained 
for some range of exponents. Since the field dependence 
of z g is very weak, and since the error bars of z g over- 
lap over some range for all the fields studied here, the 
average value z g = 3.8 ± 0.2 is accepted as the result of 
this analysis. Note that for a given z g , the corresponding 
value of Tg can be read off immediately from Fig. 1, with 
the small error of ±0.02 K. However the total error in T g , 
corresponding to the ambiguity in defining z g , is found 
to be 0.4 K for the fields used here. 

This crossing point method for identifying z g and T g 
optimizes scaling near T — T g , and involves only the 
data taken near T g . In the analysis developed above, 
the effects of the normal state conductivity have been ig- 
nored, assuming the fluctuation contributions to be much 
larger in comparison. To justify this assertion, the total 
non-linear conductivity is compared to the normal state 
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contribution at T = T g , which is estimated roughly by 
extrapolating high temperature conductivity data. For 
the current densities used here, the total conductivity is 
at least 3 orders of magnitude larger than the normal 
state conductivity. 



B. Determination of v a 

Using the results of the previous section for T g and 
z g , the static scaling exponent v g can now be obtained 
through a full scaling of the conductivity data accord- 
ing to Eq. (|l|). However, it is beneficial to form an 
extension of the crossing point method of the previous 
section, which allows a value of v g to be identified for a 
specific (correlated) pair of parameters (z g , T g ). This 
technique has the advantage of once again optimizing 
scaling near T g . To augment the crossing point method, 
a second (and more conventional) method is also used 
to determine v g . This second technique complements the 
crossing point method by optimizing scaling in the ohmic 
(high-temperature) part of the scaling region. Of the two 
different methods, the first has greater theoretical justi- 
fication, since it relies only on data in the vicinity of T g . 
However, the two methods yield nearly identical results, 
with similar error levels. This supports the fact that the 
glass critical regime is large and robust .□ We now elabo- 
rate on the two methods. 

Crossing point method, using data near T = T g : 
The method of the previous section can be extended to 
obtain terms involving v g . Since the conductivity scaling 
function F[x] of Eq. (jl]) is analytic at x — 0, a Taylor 
series expansion may be performed about this point. For 
i ~ 0, Eq. ([!]) becomes 

CT ( j/ T )(2 3 -i)/2 = F [ ] + tiTg/jy/^F'lO] + ■■■ 

= F[0] + mt + ■ ■ ■ . (2) 

Truncating Eq. (||) at its zeroth order term, F[0], gives 
the z g crossing point method used above. To determine 
v g , the first order term is also retained. The coefficient 
m = {Tg/jf/^oF'lQ] of the linear term in Eq. (g) is de- 
termined for different current densities J at T — T g by 
observing the slopes of the different curves in Fig. 1 at the 
crossing point. The result found above, z g = 3.8, is used 
on the left hand side of Eq. (g) . Since v g should be inde- 
pendent of the current, the values obtained in this way 
are averaged over the different currents to obtain error 
bars. The results of the crossing point method for v g are 
shown in Fig. 1 (inset) for each field. No obvious trend 
of v g with the field is apparent in this figure (and none is 
expected). It is therefore meaningful to average v g over 
field variations, giving a final result of v g = 1.77 ± 0.03. 
Keeping significant figures consistent with z g , we quote 
v g = 1.8 ±0.1. 

Slope method, using data in the ohmic regime: 
This makes use of the ohmic resistivity data in the limit 
of J — > and T > T g . Typically, data are obtained here 



using current densities in the range 20 to 5 x 10 2 A cm -2 . 
In the ohmic regime, p is independent of J, so Eq. (|l|) 
reduces to p ~ t Vl >( Zg ~ x > . This power law behavior can 
be observed in Fig. 2. The slope of the linear (power 
law) portion of the curve can then be used to determine 
v g (zg — 1), and subsequently v g , using z g = 3.8. The 
apparent power law behavior of p, observed over two 
decades in Fig. 2, indicates that data used to obtain v g 
in this way remain within the critical scaling region. The 
resulting v g values are shown for different fields in the 
inset of Fig. 2. The error bars for each data point reflect 
the ambiguity in identifying the limits of the ohmic scal- 
ing region. Similar to the inset of Fig. 1, no obvious trend 
of Vg with H is observed, nor is there any correlation be- 
tween the field variations in the insets of Figs. 1 and 2. 
This provides additional support for identifying v g as a 
field average: v g = 1.79 ± 0.04, a result nearly identical 
to the crossing point method. Again, keeping significant 
figures consistent with z g , we quote v g = 1.8 ± 0.1. This 
value is used in the following analysis. 



IV. DISCUSSION OF UNIVERSALITY 



The absence of field dependence in the exponents v g 
and z g , demonstrated in Section III, is consistent with 
the vortex glass theory, which states that the critical be- 
havior of the linear or nonlinear conductivity a near the 
transition line T g (H) should belong to a single univer- 
sality class, for some range of H > O.til The vortex glass 
excitations are thought to be formed within a medium 
of field-induced vortices, whose excursions from perfect 
rigid rod geometry, parallel to H, are affected by entropy, 
applied currents, and the pinning effects of quenched 
disorder. The density of field vortices varies with H , 
while quenched disorder does not. For a particular sam- 
ple then, vortex matter at different fixed fields forms 
unique vortex systems. The behavior and response of 
these unique systems to external probes may vary greatly 
with the field. Indeed, it has been noted that changes in 
the field may be responsible for phenomena such as di- 
mensional crossover, causing vortex matter at low and 
high fields to behave very diffcrently.E2l However, many 
of these phenomena are of mean-field origin, and should 
not be observed in measurements which predominantly 
reflect a critically diverging correlation length. It is not 
surprising then that vortex glass universality, observed 
here for all studied fields, H > 0, along T g (H) (and in 
Rcf. 15 up to much higher fields) appears to hold in spite 
of such variations with field. 

In order to observe universal glass behavior of both 
the critical exponents and the scaling function, it is nec- 
essary to introduce a new notation to account for field 
dependence in the vortex glass description. This field 
dependence is not specified by the vortex glass theory; it 
is not known to be universal, a priori, and may vary from 
sample to sample. The conductivity scaling of Eq. (g) can 



3 



then be rewritten as 



1*1 



-V g {Zg-l) 



J T 

- — 1*1 
T Jo 



-2v Q 



(3) 



Here, we have introduced the field dependent resistivity 
scale po(H) and the current scale Jo (if). The transition 
line T g (H) is also field dependent. The new functions 
G±[x], corresponding to T > T g and T < T g , and their 
arguments, are dimensionless. 

To further develop a universal description, it is neces- 
sary to identify the scales po(H) and Jq{H) empirically. 
Methods for obtaining characteristic resistivity and cur- 
rent scales exist in the literature ,fl however we believe the 
following definitions to be the most objective: 



Po = pt 



-Ug(Zg-l) 



T>T g ,J-*0 

Jo = J(p/po) 2/{1 - z ^ 



T=T„ 



(4) 
(5) 



where p and J refer here to experimentally determined 
quantities. In this definition, po is the ohmic resistiv- 
ity scale, while Jo is the current required to produce a 
resistance p = po along the critical isotherm. 

Using these definitions, the nonlinear resistivity curves, 
scaled at different fields according to Eq. , collapse as 
shown in Fig. 3. The choice of normalization in Eqs. (^) 
and causes the dimensionless point (1, 1) to be located 
as shown in the figure. We emphasize that this collapse 
of J-E scaling curves, including different fields, demon- 
strates the universality of the vortex glass transition. By 
varying v gi z g , and T g (H) within their error bars, as 
found above, it is not possible to obtain any discernible 
improvement in the scaling beyond Fig. 3. Additionally, 
for z g < 3.8, the scaling deteriorates markedly. 

The obtained field dependences of po and Jo are shown 
in Fig. 4. As mentioned above, neither of these quantities 
is specified by the vortex glass theory. However, they can 
be determined in a multicritical theory such as 3D XY, 
which is thought to apply to topological fluctuations at 
low fields Jj'EJEJ The following power law dependencies pf 
Po{H) and Jo (if) are derived from the 3D XY theory x3 



J (H) = JH 1 -^ U "\ 



-1)1/21/* 



(G) 
(7) 



test is propded by Eq. (^), using the known value of 
v xy = 0.67.E3 The obtained value of v g = 1.60±0.03 ~ 1.6 
provides a satisfactory check on the universality analysis. 
Additionally, z xy may be calculated from Eq. (0), using 
the slope obtained in Fig. 4, together with v xy = 0.67, 



1.8, and z g = 3.8. The result is z x 



2.3 ± 0.2. 



[Note that the error sources in Eq. (g) are correlated.] 
This result is closer to that given by Booth et a/.Ej 
than Moloni et al&B We speculate that the difference be- 
tween the present results and those of Ref. |l5| may reflect 
the difficulty in treating background contributions to a, 
which impedes the analysis of Ref. |lj, but not the present 
one. 



V. CONCLUSIONS 



In this work, the vortex glass scaling parameters v g = 
1.8 ± 0.1, z g = 3.8 ± 0.2, and T g (H) have been obtained 
through a systematic scaling technique, involving the vor- 
tex glass crossing point at T = T g . Similar to other de- 
oxygenated YBa2Cu307_5 films, the vortex glass scaling 
region was found to be wide, on the order of 10 K at 
H = 7 T. The absence of field dependence in the expo- 
nents v g and z g provides a demonstration of universal 
behavior, since different fields correspond to unique vor- 
tex systems. A second (and stronger) demonstration of 
universality is given by the collapse of scaled J-E curves 
obtained at different fields in the vortex glass critical re- 
gion, as shown in Fig. 3. Because universality is an es- 
sential component of critical phenomena, these checks on 
both the glass exponents and the scaling functions pro- 
vide strong evidence for a critical phase transition line 
along T g (H) for H > 0. While additional evidence has 
not been presented here to demonstrate that this transi- 
tion involves explicitly the broken symmetry associated 
with a vortex glass, it should be noted that the scal- 
ing functions in Fig. 3, below and above T g (H), exhibit 
the expected glassy characteristics at low currents: below 
Tg(H) the non- linear resistivity decays exponentially to 
0, while above T g (H) the resistivity decays exponentially 
to its normalized value of 1. 

The values reported here for the dynamic scaling expo- 
nent, z gi and the static scaling exponent_j/ g , fall within 
the range of previously reported values JIj although our 



where v xy and z xy are 3D XY scaling exponents, and p result-.z 9 = 3.8 ± 0.2 falls barely within the expected 



and J are unknown constants which are sample-, but not 
field-dependent . 

In Fig. 4, po(if)r a pd Jo (if) are described well by power 
laws for low fieldsO It is important to note that, in con- 
trast to other 3D XY analyses performed in the liter- 
ature, the observance of this power law does not hinge 
on an appropriate choice of T c as a fitting parameter. 
In the fits shown, only data are used which fall within, 
the expected 3D XY scaling region for this samplcE£l 
(H < 8 T). We find corresponding slopes of -3.1 ± 0.2 
(po) and -1.39 ±0.04 ~ -1.4 (J ). A self-consistency 



ranged z g > 4. The significant variation of reported z g 
values in the literature could imply that conventional 
methods for determining z g may involve systematic er- 
rors in the identification of the critical isotherm (and thus 
T g ), a problem over which we believe the techniques ap- 
plied here provide greater control. Alternatively, the field 
variations of z g could result from the glass correlation 
length becoming jSmaller than the distance between field- 
induced vortices. Q'tll Finally, z g variations could imply an 
absence of dynamic universality. These possibilities are 
not investigated in further detail here. However, to study 
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this issue in more depth, it will be important to apply the 
methods developed in this work to other high-T c samples. 
In each case, universality must be demonstrated for both 
the scaling exponents and the functions. 

Finally, through Eqs. (^) and (]?]), evidence has been 
presented which, together with analyses such as Ref. |l^, 
provides strong support for 3D XY multicritical scal- 
ing. The present analysis involves scaling of the nonlinear 
resistivi ty , an d therefore complements previous 3D XY 
studies,OlIZE3 which apply only to the linear resistiv- 
ity. The focus of the present work has been the vor- 
tex glass critical regime, for which Eqs. (§),(§), and (0) 
augment previous analyses by explaining the field depen- 
dence of vortex glass scaling, in addition to the temper- 
ature dependence. Through these equations it is found 



that z xy ~ 2.3, in agreement with Ref. 23; the present 



result does not suffer from background subtraction diffi- 
culties. Eqs. (|J) and (^) apply only to the 3D XY scaling 
region, thus allowing the range of 3D XY scaling to be 
determined. For the sample studied here, the 3D XY 
critical regime extends to 8 T along the line T = T g (H). 
If the dynamics of the glass transition are ultimately de- 
termined to be non-universal, it will be important to as- 
certain whether the self-consistency of 3D XY multicrit- 
ical scaling, demonstrated here, still persists. 
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FIG. 1. Crossing point method at H = 1 T. Differ- 
ent curves correspond to different current densities. The 
curves all cross at a single point, T — 65.2 K, provided 
z g is chosen correctly, thus identifying T g . The z g val- 
ues obtained in this way are shown in the inset for each 
applied field. Additionally, by determining the slopes of 
the curves at T g , the exponent v g can be obtained, as 
described in the text. These v g values are shown for each 
applied field in the inset. 
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FIG. 2. Slope method for H = 1 T. The static scaling 
exponent v g is calculated from the power-law dependence 
of the ohmic p(T) data. Values of T g and z g — 3.8 used 
in this fit were obtained from the crossing point method 
of Fig. 1. Resulting values of v g are shown in the inset 
for each applied field. 

FIG. 3. Scaled p(J, T) data from three different fields, (p 
and J represent, respectively, the left hand side of Eq. (3), 
and the argument (x) of the scaling function G± [x] . Di- 
mensionless, arbitrary units are used.) The quantities po 
and Jo of Eq. (3) have been found for each field so that 
the three curves collapse onto a single, universal curve. 
The scaling exponents z g = 3.8 and v g = 1.8 are used for 
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all fields. Although the collapse is excellent for all five 
fields used here, only three are shown for clarity. The 
dashed lines represent asymptotic behaviors, and their 
crossing identifies the normalized point (1, 1); see Eqs. (4) 
and (5). 



FIG. 4. Field dependence of the quantities po and 
Jo- The dashed lines represent best fits, involving data 
points only from the expected 3D XY scaling regionE-2 
(H < 8 T). 
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